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ABSTRACT 

With the advent of the Event Horizon Telescope (EHT), a millimeter/sub-millimeter very-long baseline interfer- 
ometer (VLBI), it has become possible to image a handful of black holes with sub-horizon resolutions. However, 
these images do not translate into microarcsecond absolute positions due to the lack of absolute phase information 
when an external phase reference is not used. Due to the short atmospheric coherence time at these wavelengths, 
nodding between the source and phase reference is impractical. However, here we suggest an alternative scheme 
which makes use of the fact that many of the VLBI stations within the EHT are arrays in their own right. With 
this we show that it should be possible to absolutely position the supermassive black holes at the centers of the 
Milky Way (Sgr A*) and M87 relative to nearby objects with precisions of roughly 1 /ias. This is sufficient to 
detect the perturbations to Sgr A*'s position resulting from interactions with the stars and stellar-mass black holes 
in the Galactic cusp on year timescales, and severely constrain the astrophysically relevant parameter space for 
an orbiting intermediate mass black hole, implicated in some mechanisms for producing the young massive stars 
in the Galactic center. For M87, it allows the registering of millimeter images, in which the black hole may be 
identified by its silhouette against nearby emission, and existing larger scale radio images, eliminating present 
ambiguities in the nature of the radio core and inclination, opening angle, and source of the radio jet. 
Subject headings: Black hole physics — Galaxy: center — Techniques: interferometric — Submillimeter: general 
— Astrometry — Proper motions 



1. INTRODUCTION 

The Event Horizon Telescope (EHT), a proposed Earth-sized 
array of existing millimeter and sub-millimeter observatories, 
promises to provide microarcsecond imaging resolutions via 
very-long baseline interferometry (VLBI) (Doeleman et al. 
2009a). This is sufficient to resolve the horizons of a hand- 
ful of supermassive black holes, including Sagittarius A* (Sgr 
A*) and M87. Already, a mm- VLBI array consisting of three 
stations has probed sub-horizon scale structure in Sgr A* at 
1.3mm wavelengths (Doeleman et al. 2008; Fish et al. 2011). 
These have been used to test fundamental physics, providing 
direct evidence for the existence of an event horizon surround- 
ing the supermassive black hole at the center of the Milky Way 
(Broderick et al. 2009b), as well as study the accretion flow 
surround Sgr A*, giving an estimate for the black hole spin 
(Broderick et al. 2009a; Huang et al. 2009; Moscibrodzka et al. 
2009; Dexter et al. 2010; Broderick et al. 2010). 

The existing observations have made use of the James Clerk 
Maxwell Telescope (JCMT) and the Sub-Millimeter Array 
(SMA), located atop Mauna Kea in Hawaii, the Combined Ar- 
ray for Research in Millimeter-wave Astronomy (CARMA) 
in Cedar Flat, California, and the Arizona Radio Observa- 
tory Sub-Millimeter Telescope (SMT) on Mount Graham, Ari- 
zona, giving a maximum baseline of 4500km, oriented pre- 
dominantly east-west. Potential future stations include sites in 
Chile (Atacama Pathfinder Experiment, Atacama Submillime- 
ter Telescope and Atacama Large Millimeter Array; APEX, 
ASTE, and ALMA, respectively), Mexico (Large Millime- 
ter Telescope; LMT), the South Pole (South Pole Telescope; 
SPT), and the IRAM telescopes in Spain (Pico Veleta; PV) 
and France (Plateau de Bure; PdB). Among these the longest 
baselines are ~ 1.2 x 10 4 km, corresponding to a synthesized 
beam of 30 /ias at 230GHz (1.3mm) and 20 ^as at 345GHz 



(0.87 mm). 

In addition to including new baselines, improvements in sen- 
sitivity are actively being pursued, both by pushing towards 
larger bandwidths and collecting areas. Using the JCMT, SMT 
and a single 10m antenna at CARMA, Fish et al. (2011) was 
able to measure 0.1 Jy visibilities on long baselines. Efforts 
already underway to increase the bandwidth to 4 GHz will im- 
prove that by a factor of two. Phasing together the multiple 
telescopes on Mauna Kea (SMA, JCMT & CSO; Weintroub 
2008) and many antenna at CARMA will improve this further, 
increasing sensitivity by a factor of 3 at some sites, making it 
possible to detect sources considerably dimmer than Sgr A* 
and M87. This will make it possible to contemplate doing as- 
trometry on microarcsecond scales using the EHT. 

Since the advent of self-calibration, phased-reference VLBI 
observations have only been necessary in imaging observations 
for low-brightness sources. Nevertheless, because the use of 
closure phases discards an over-all phase calibration, phase- 
referenced observations provide the unique capability to do 
relative astrometry between distinct sources, e.g., positioning 
a target relative to distant background quasars. In the case of 
the Galactic Center, 7 mm- VLBI observations using the Very- 
Long Baseline Array have already paid substantial scientific 
dividends: placing tight constraints upon the motion of Sgr A* 
(Reid et al. 1999a; Reid & Brunthaler 2004; Reid et al. 2008), 
providing independent parallax measurements of our distance 
from the Galactic center (Reid et al. 2009), and registering 
the radio and infrared reference frames of the central 0.1 pc 
(Menten et al. 1997; Reid et al. 2003, 2007). All of these make 
specific use of background quasars to locate sources within the 
Galactic center to sub-milliarcsecond precisions. 

Moving towards shorter wavelengths is challenging for a 
number of reasons, including the typically rapidly declin- 
ing source flux densities and shortened atmospheric coher- 
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ence time (10s) and angle (isoplanatic angle, ~ 3°, see Ap- 
pendix A). Nevertheless, there are strong scientific motivations 
for attempting to do so. Chief among these is that at mm- 
wavelengths the detection of the black hole silhouette unam- 
biguously identifies and locates the supermassive black holes 
in Sgr A* and M87. Registering the mm and cm-VLBI ref- 
erence frames then provides a means to determine where the 
black hole is relative to larger morphological features. Mi- 
croarcsecond astrometry also allows the direct detection of the 
jitter induced by the interaction of the supermassive black holes 
with the retinue of compact objects in their vicinity, providing 
a probe of the dark component of the stellar cusp surrounding 
the central object. Finally, absolute astrometry makes it possi- 
ble to study short-timescale morphological changes in Sgr A* 
and M87. 

Here we describe some of the the scientific motivations for 
phase-referenced mm-VLBI in more detail, giving estimates of 
characteristic time scales and precisions needed to detect the 
relevant phenomena. Since the particular questions that are ac- 
cessible are different in Sgr A* and M87, we will address each 
separately. Section 2 describes a possible scheme by which 
mm-VLBI may be accomplished, including some of the criti- 
cal limitations that the short wavelength imposes. Sections 3 
and 4 provide the motivations as well as candidate phase refer- 
ences that meet the requirements set out in Section 2. Finally, 
concluding remarks may be found in Section 5. 

2. PHASE REFERENCED MILLIMETER- VLBI 

With the advent of self-calibration, interferometric radio ob- 
servations became capable of reconstructing the vast majority 
of complex visibility phase information in a given observation, 
correcting for the nearly all of the station-dependent errors due 
to propagation through the atmosphere. These methods are, 
however, insensitive to an overall phase shift of the complex 
visibilities, including to a translational uncertainty in the im- 
age. This precludes VLBI observations from measuring the 
absolute position of radio features at the incredible precision 
that it can image them. 

This situation may be ameliorated with the use of a phase 
reference, which provides a fixed point against which to mea- 
sure the relative location of the source. For this purpose, back- 
ground quasars are frequently used, due to their small intrinsic 
angular size, high brightness and exceedingly small transverse 
velocities. In practice, the small angular size and fixed loca- 
tion of the phase reference provides an easily modeled system 
with which to reconstruct with high accuracy the phase errors 
at each VLBI station, fulfilling an analogous role as guide stars 
in optical and infrared adaptive optics systems 1 . 

Any potential phase reference must be sufficiently close 
to the target that it suffers strongly correlated phase distor- 
tions, i.e., it must lie within the isoplanatic patch. At mm- 
wavelengths this depends upon the particular site, typically 
requiring separations < 2°-3°, though we will discuss this 
in more detail below (see also Appendix A). Similarly, the 
phase reference must be detectable on timescales compara- 
ble to the atmospheric coherence time, roughly 10 s at mm- 
wavelengths. In practice, these impose severe limitations 
on phase-referenced mm-VLBI. While a handful of potential 
phase references may exist within 2° of the target source, none 

1 In this context, when we refer to "phase-referenced mm-VLBI" we are not 
describing a scheme for performing mm-VLBI observations on dim sources. 
Rather, we are discussing a means to assess the relative positions of sources 
that are sufficiently bright to be detected independently. 



Table 1 

EHT Baseline Sensitivities for Multi-Antenna Stations 



Realization Baseline f23o(mJy) ^345(mJy) '2300^-' ^ilanQy) 



EHT-I 


H J -C' 


41 


63 


1.3 


1.7 




rf-A 1 


19 


34 


0.61 


0.91 




C'-A 1 


27 


45 


0.88 


1.2 




P'-A 1 


17 


36 


0.56 


0.94 


EHT-C 


H a -C 8 


10 


16 


0.33 


0.42 




H a_ A 10 


3.9 


7.1 


0.13 


0.19 




C 7 -A 10 


3.3 


5.5 


0.11 


0.15 




p6_ A 10 


2.2 


4.6 


0.072 


0.12 



lie within the beam widths of individual EHT stations, which 
range from 6" (LMT) to 54" (SMA). As a consequence, ob- 
servations of the target source and phase reference with a sin- 
gle dish requires repointing the telescopes. However, this must 
be accomplished well within the atmospheric coherence time, 
leaving sufficient time to detect both the source and the phase 
reference, something that is presently not possible. 

However, the EHT is unique in that many of the VLBI sta- 
tions are interferometers in their own right, i.e., contain many 
independent antennas. As a consequence, at these stations it 
is possible to observe both the target source and phase refer- 
ence simultaneously using a subset of the available antennas 
for each (the Paired Antennas approach described in Chapter 
13.2 of Thompson et al. 2001). The remaining instrumental 
components of the complex gains change slowly in compari- 
son to the atmospheric propagation effects, and thus may be 
measured by periodically switching antennas from one sub- 
set to the other. This is identical to the approach used by 
Counselman et al. (1974) to detect the gravitational deflection 
of quasars via VLBI, and similar to the principles underlying 
VERA (Kawaguchi et al. 2000; Honma et al. 2000). 

Stations consisting of arrays are some of the most sensitive in 
the EHT, providing the most permissive flux limits upon poten- 
tial phase references. Potential stations for this purpose include 
Hawaii, CARMA, ALMA and Plateau de Bure. The limiting 
sensitivities at 230 GHz and 345 GHz on baselines that include 
these stations at both ends, assuming a bandwidth of 4 GHz and 
10 s integration times, are listed in Table 1. We include sensi- 
tivity estimates for two possible realization: 

EHT-I: Using only the JCMT in Hawaii (H J ), a single 10.4m 
antenna at CARMA (C 1 ), 12 m antenna at ALMA (A 1 ), 
and a single 15 m antenna at Plateau de Bure (P 1 ). 

EHT-C: Using the phased Hawaii station (including the SMA, 
JCMT and CSO, H a ), 8 CARMA antennas (6 10.4m and 
2 6.1m antennas, C 8 ), 10 12m ALMA antennas (A 10 ) 
and all 6 15 m Plateau de Bure antennas (P 6 ). 

The first of these indicate what is possible with limited num- 
ber of antenna, similar to what is possible for initial version of 
the EHT, available in the near future. The second corresponds 
to the maximum achievable sensitivities using the complete 
EHT though employing only existing facilities. To estimate 
these flux limits we used the system-equivalent flux densities 
reported in Doeleman et al. (2009b) 

In Appendix A we estimate the atmospheric limits upon the 
accuracy with which angular separation may be measured. In 
principle, this is also limited by source structure and variability, 
though imaging with the EHT should ameliorate these (indeed, 
as described in Section 3.1.4 the characterization of source 
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variability is one of the motivations for high-precision astrom- 
etry), and instrumental phase noise. Atmospheric limitations 
come in two forms: that due to the small-scale, rapidly vary- 
ing structures that limit the coherence timescale, and that due 
to large-scale, slowly evolving features that induce anomalous 
phase delays resulting from the slighly different zenith angle, z, 
of the target and reference sources. As a result, the mitigation 
strategies for these differ substantially. Generally, we find that 
sub-beam precision is almost certainly achievable, with /ias ac- 
curacy likely possible. Note, however, that this accuracy is not 
always required. For example, the precision with which it is 
possible to register mm-VLBI and existing cm-VLBI images 
is necessarily limited by the resolution of the latter, which is 
roughly 0.1 mas. 

As described in Appendix A, the short timescale of the rapid 
component implies that hundreds of realizations of the atmo- 
spheric turbulence are encountered in a typical observing night 
(roughly 2hr of on-source integration time). For phase refer- 
ences within 1 .4° the atmospheric conditions will permit beam- 
scale (30 /ias) astrometry on 10 s intervals at 230 GHz. Over 
10 min, the integration times used in Doeleman et al. (2009b) 
and Fish et al. (201 1), this can be reduced to 4 /ias scales. In a 
2 hr period this may be reduced to 1 /ias. Thus, in practice, this 
is unlikely to be the dominant contribution to the astrometric 
uncertainty. 

At 7 mm, the slowly varying, large-scale component domi- 
nates the astrometric uncertanties (Reid et al. 1999a), produc- 
ing typical unmodeled time delays of 0.1ns for sources sepa- 
rated by 1 in azimuth. This corresponds to an uncertainty on 
the order of 0.1 mas, i.e., three beam widths at 1.3 mm. How- 
ever, there are a variety of schemes that may be employed to 
reduce this significantly. All of these seek to do so by exploit- 
ing the large spatial and temporal scales of the atmospheric 
fluctuations responsible for this component of the astrometric 
uncertainty. 

The first is to implement better atmospheric models. The 
atmospheric model of Niell (1996) already accounts for the 
seasonal variations, generating daily estimates for the zenith 
phase delay, and thus does model the very-long timescale fluc- 
tuations. Improvements over this model come in a variety 
of forms, including the structure of the atmosphere assumed 
(Boehm et al. 2006a) and the inclusion of real-time weather 
modeling to produce accurate maps on 6 hr timescales (Boehm 
et al. 2006b). These promise to reduce the size of the unmod- 
eled phase perturbations by factors of a few. 

The second is to periodically measure the group delays di- 
rectly. This can be done using sets of known bright quasars, 
which are assumed not to move, so-called geodetic blocks, and 
has already been done in astrometric studies of the Galactic 
Center (Reid & Brunthaler 2004; Reid et al. 2009). Alter- 
natively, this can be done using Global Positioning Systems 
(GPS) and Global Navigation Satellite Systems (GNSS) (see, 
e.g., Byun & Bar-Sever 2009). Already, GPS/GNSS is used to 
produce estimates for the zenith path delay with a time resolu- 
tion of 5 min and precision of roughly 3 mm at sites comparable 
to those that will be employed by the EHT This alone would 
be sufficient to improve the astrometric uncertainty due to the 
slow component by an order of magnitude, giving an accuracy 
of roughly 10 /ias. 

If three or more viable references are located in the vicinity 
of the target 2 , a third strategy is to calibrate the unknown large- 

2 Note that because the atmospheric fluctuations that generate the slowly- 
varying phase perturbations are necessarily large, additional reference sources 



scale phase delays directly. This can be done by asserting that 
the true positions of the references are fixed, and then using 
a procedure similar to that described in Appendix B to esti- 
mate the site specific phase delays by comparing the observed 
and expected reference separations (see also Fomalont 2005). 
This assumes that the references are not variable themselves, 
though even at /ias resolutions distant quasars are unlikely to 
exhibit structural variability on hour timescales. In principle, 
as shown in the appendix, the accuracy with which this can be 
done is comparable to the precision limit associated with the 
short-timescale variability, roughly 1 /ias. 

Finally, the slow variations can be averaged down by com- 
bining observations over many days. The upper limit upon the 
time over which the observations can be combined is dictated 
by the time-scale of variations that are of interest. Thus, e.g., 
fluctuations in the position of Sgr A* on year-long timescales 
allows averaging over ~ 3 months, corresponding to ~ 10 2 in- 
dependent realizations of the large-scale fluctuations, and thus 
a reduction in the uncertainty due to the slow-component of an 
additional order of magnitude. 

In practice a combination of these mitigation strategies are 
likely to be required. Nevertheless, there is clearly a good 
reason to believe that /ias accuracies are achievable over long 
times via improved atmospheric modeling, direct phase delay 
measurement and time averaging. Furthermore, it may be pos- 
sible to reach these accuracies over much shorter times using 
multiple phase references, though this is more speculative. 

In the following sections we discuss some of the science 
made possible by phase-referenced mm-VLBI observations, 
and suggest possible phase references to be used for these pur- 
poses. Since many of the scientific motivations are different for 
Sgr A* and M87, the observing strategies also differ. There- 
fore, we address these sources separately. 

3. SAGITTARIUS A* 

Associated with the supermassive black hole located at the 
center of the Milky Way, Sgr A* is the primary target for 
the EHT. Presently, the best estimates of the mass (Af) and 
distance (D) of Sgr A* come from the observations of orbit- 
ing stars (the so-called "S-stars"). These have yielded M = 
4.3 ± 0.5 x 1O 6 M and D = 8.3 ± 0.4kpc, respectively, where 
both include the systematic uncertainties (Ghez et al. 2008; 
Gillessen et al. 2009b, a). This results in an angular size for 
the apparent horizon of 53 ± 2 /ias 3 , the largest for any known 
black hole. The mass of Sgr A* is necessarily confined to 
within the periapse of nearby stars, giving a maximum radius 
of roughly 10 2 AU ~ 3 x 10 3 GM/c 2 . 

3.1. Scientific Objectives 

A variety of astrometric studies of Sgr A* have already been 
proposed and carried out. As a consequence of the large extinc- 
tion towards the Galactic center, these necessarily are confined 

are not subject to the same constraints that limit the reference-target separa- 
tion. Since the distortion in the separation grows linearly with zenith separa- 
tion, larger separations between reference sources result in larger relative phase 
perturbations due the large-scale fluctuations. However, the uncertainty associ- 
ated with the rapid variations also grow nearly linearly, and thus the fractional 
precision with which the slowly-varying phase perturbations can be measured 
grows only slowly. Nevertheless, the number of bright references grows as sep- 
aration squared. The separations are fundamentally limited by the region over 
which the large-scale atmospheric fluctuations can be accurately modeled. 

3 Because the current mass and distance measurements are strongly corre- 
lated, with mass scaling roughly as M oc D 1 ' 8 (Ghez et al. 2008), the angular 
size of Sgr A* is constrianed much more strongly than M or D alone. 



4 



BRODERICK, LOEB & REID 



to the radio and near-infrared (NIR). Nevertheless, a wide va- 
riety of objectives have been and are being pursued, including 
placing independent constraints upon the mass of the central 
supermassive black hole (Reid et al. 1999b; Reid & Brunthaler 
2004, 2005) and detecting structural variations during radio 
and NIR flares (Reid et al. 2008; Bartko et al. 2009). A generic 
feature of these is the need for high astrometric precision. Here 
we describe how these efforts may be complemented or im- 
proved upon using the /ias precision afforded by mm- VLBI, as 
well as suggest additional scientific motivations for developing 
a mm- VLBI astrometry capability. 

3.1.1. Constraints upon the Proper Motion of Sgr A * 

The current hierarchical paradigm for black hole formation 
and active galactic nuclei (AGN) fueling is fundamentally dy- 
namic, imparting large velocities to the central supermassive 
black holes. However, on timescales short in comparison to the 
time between galactic mergers this motion is dissipated by col- 
lective interactions between the supermassive black hole and 
the stars in the stellar bulge, i.e., dynamical friction (Chatter- 
jee et al. 2002). Thus we expect the black hole at the center of 
the Milky Way to be nearly at rest, located at the bottom of the 
Galactic potential well. However, because of stochastic inter- 
actions with nearby stars, it will not be completely stationary 
and, if the surrounding stellar core is well described by a Plum- 
mer sphere, will execute random motions (perpendicular to the 
Galactic plane, say) with variance 



2 GM tot 
9 a ~M 



(1) 



where M tot is the total mass of the nuclear stellar cluster, m 4 
is the characteristic mass of a star and a is the core radius of 
the Plummer sphere (Chatterjee et al. 2002; Reid & Brunthaler 
2004). Thus, given an estimate of the stellar content within 
some radius, R, it is possible to constrain the the mass of Sgr 
A* given an upper limit upon its motion: 



M > 



M R 



1 + (9 /2) { (vi)aR 3 / \G (R 2 + a 2 ) 3/2 «*] } 



(2) 



(Chatterjee et al. 2002; Reid & Brunthaler 2004). Note that this 
constraint is independent of those obtained from measurements 
of the orbits of stars immediately surrounding Sgr A*. 

High-precision astrometric studies of Sgr A* have already 
been performed at 7mm using the Very Long Baseline Ar- 
ray (VLBA), constraining (v^) < 2kms _1 , corresponding to 
M > 2 x 1O 6 M0. They have also measured the apparent veloc- 
ity parallel to the Galactic plane, giving vii ~ 250± 16kms _1 
(assuming a distance of 8.3 kpc to Sgr A*), which is presum- 
ably dominated by the contribution from the Sun's orbit in the 
Galaxy. The VLBA observations have a nominal astrometric 
precision of 10 2 /ias, and extend over a period of more than a 
decade. During this time they see no evidence for fluctuations, 
implying that either the motion is secular or has a period sig- 
nificantly exceeding ~ 20 yr. Note that the natural period for 
oscillations about the center of the average collective Galactic 
potential is dependent upon the radius (/?j n f) of the region of in- 

fluence, and roughly 7 .4 x 1 3 (Rm /0 . 3 pc) yr, much larger 
than any observational period of interest. 

Empirical estimates for the magnitude of the secular motion 
of Sgr A* are presently limited by the accuracy of individual 
radio observations. In principle, the uncertainty in the perpen- 
dicular component of the peculiar motion of the Sun provides 



a fundamental limit upon efforts to measure the secular motion 
of Sgr A*. This is presently estimated to be 7.25 ± 0.37kms _1 
relative to the local standard of rest (Schonrich et al. 2010). 
An additional uncertainty may arise from Galactic precession, 
which can be as high as 0.6 km s" 1 in the extreme case in which 
the Milky Way "flips" in a Hubble time. These corresponds 
to an angular deviation of 9.4 /ias yr" 1 and 15/iasyr" 1 , respec- 
tively, both much smaller than the limits imposed by the ac- 
curacy of existing radio observations. Therefore, considerable 
improvement in the astrometric accuracy is possible in princi- 
ple. 

3.1.2. Detection of the Evolved Component of the Cusp 

In addition to a secular component due the Sun's orbit and 
a potential displacement from the center of the Galactic poten- 
tial, the velocity of the central supermassive black hole has a 
fluctuating component due to random forces applied by fluc- 
tuations in the surrounding retinue of stars and their compact 
remnants. In contrast to the secular motion, detecting a fluctu- 
ating component is made much simpler due to the variations in 
the velocity, which clearly differentiates it from contributions 
due to the solar motion. With /ias astrometric precision, detect- 
ing this fluctuating component on scales of a few years should 
be possible. 

It is from this fluctuating component that the mass estimates 
described in the preceding section is ultimately derived. How- 
ever, it is also possible to probe the spatial and mass distribu- 
tion of the stars and their remnants in the immediate vicinity of 
Sgr A*. Understanding the properties of these objects is of in- 
trinsic interest; how stars form and evolve in the Galactic center 
is presently poorly understood, characterized by the so-called 
"paradox of youth" (Ghez et al. 2003). However, it is also im- 
portant as an exemplar. Detecting the stellar component would 
provide a direct means with which to estimate the rates of stel- 
lar disruptions and the subsequent contribution to the growth 
of AGN. 

Furthermore, orbiting populations of stellar-mass black holes 
are expected to play a pivotal role in present and future grav- 
itational wave searches. Stellar-mass black hole binary coa- 
lescence rates for the Laser Interferometer Gravitational Wave 
Observatory (LIGO) are strongly affected by the concentration 
of the stellar-mass black holes in galactic nuclei (O'Leary et al. 
2009). Similarly, the event rate of extreme mass-ratio inspirals 
(EMRIs), one of the primary science drivers for the Laser Inter- 
ferometer Space Antenna (LISA) (see, e.g. Miller et al. 2009; 
Schutz et al. 2009), is strongly dependent upon the density and 
mass function of stellar-mass black holes near the supermas- 
sive black hole (Sigurdsson & Rees 1997; Gair et al. 2004; 
Hopman & Alexander 2005, 2006; Gair et al. 2010; Amaro- 
Seoane & Preto 2010). We note that the environment of Sgr A* 
is particularly relevant to LISA, since it is comparable in mass 
to the systems that will dominate the LISA frequency band. In 
both cases, empirical estimates of the black hole distribution in 
Sgr A* can substantially reduce the uncertainties in the inferred 
merger rates. 

Here we estimate the precision necessary to directly observe 
the motion induced upon Sgr A* by the nuclear stellar and 
black hole populations. To do this we assume a fiducial set 
of distributions for these objects, and consider both optimistic 
and pessimistic deviations from these. In all cases the natu- 
ral frequency of oscillation about the minimum of the Galactic 
potential and the damping timescale due to dynamical friction 
with the stellar bulge (ljq and 7 in Appendix C) are relevant 
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only at much larger timescales than the observationally moti- 
vated subset we consider (i.e., sub-century). 

The stellar mass density surrounding Sgr A* has been in- 
ferred from infrared studies, finding that the total stellar mass 
density is well approximated by a broken power law: 



p*(r)~ 1.7 ±0.8 x 10 6 



0.22 pc 



M Q pc" 



(3) 



where 7 = 1.2 for r < 0.22pc and 7 = 1.75 for r > 0.22 pc 
(Schodel et al. 2007). This is well constrained observationally 
for r > 10~ 2 pc, though there is considerable uncertainty regard- 
ing the structure of the core on smaller scales (Merritt 2010). 
In addition there is evidence for a preponderance of young stars 
in the vicinity of the supermassive black hole (Buchholz et al. 
2009; Do et al. 2009, 2010). This biases the stellar mass func- 
tion in the Galactic center towards higher masses. Neverthe- 
less, here we make the conservative choice of adopting a stellar 
mass function indicative of the Galactic disk. Specifically, we 
employ the log-normal mass function for low-mass stars de- 
termined by Covey et al. (2008), smoothly joining a Salpeter 
high-mass tail, i.e., (j)(m,r) oc to -2,3 above m ~ O.6M , and 
truncated at 1O 2 M . For our purposes here, it is enough to note 
that the average stellar mass and it's square are = 0.75 M 
and /i 2 ~ 13.7M Q 2 , respectively. Note that the Salpeter mass 
function immediately implies that /i 2 (r) is determined by the 
maximum mass and is generally much larger than m 2 . 

Less well constrained is the contribution from stellar rem- 
nants. Large numbers of stellar-mass black holes are expected 
to collect via dynamical friction in the vicinity of the central su- 
permassive black hole (O'Leary et al. 2009; Alexander & Hop- 
man 2009). A number of authors have estimated the number 
of stellar-mass black holes in the Galactic center, finding that 
roughly 2 x 10 4 are expected to exist within the inner pc (Mor- 
ris 1993; Miralda-Escude & Gould 2000; Hopman & Alexan- 
der 2006; Freitag et al. 2006; Alexander & Hopman 2009). The 
distribution of these remnants is less well understood; predic- 
tions ranging from flat cores to power-law cusps with indexes 
as extreme as -2.75 can be found in the literature (see, e.g., 
Alexander & Hopman 2009; Merritt 2010), though cusps with 
an indexes near -2 are most common. Thus we adopt a fiducial 
model for the black hole component given by 



4n 



M pc 3 , 



(4) 



with m, = IOMq and N pc is the number of black holes in the 
inner pc. 

In Appendix C we derive an approximate power spectrum 
for the displacement of the central supermassive black hole, 
Px,u>, for an arbitrary number density of objects, n(r), and a 
potentially radially dependent mass function, <j>(m,r): 



-4 



dr 



2(2tt) 3 / 2 GV n 
3 



(5) 



with ilk given by Equation (C10). In this /i is the RMS mass 
at a given radius (thus weighted towards the more massive ob- 
jects). At this point, there are some general features of note. 
First, since P x u depends upon the object mass function through 
/j, 2 n ~ fj?p/m alone, when the same mass density is present, 
perhaps unsurprisingly, more massive stars/remnants result in 
proportionally larger displacements. Second, if /x 2 is fixed, 
n(r) ex r~ v with v > 1/2, and the motions of the stars/remnants 




Figure 1. Characteristic displacements as a function of time-scale for when 
stellar mass black holes and stars are considered (solid), only the stellar com- 
ponent is considered (dot-dashed), when the stellar component is assumed to 
form a core of size 0.01 pc (long-dashed), and when all components are trun- 
cated inside of 0.01 pc (short-dashed). The region excluded by 7mm observa- 
tions with the VLBA is shown by the hatched area (Reid & Branthaler 2005). 
Along the left and right axes the displacements are shown in units of black hole 
mass and angle. Along the top axis the radius with Keplerian period equal to 
the time period is listed, providing some sense for the location of the objects 
that are dominating the fluctuations. 

is dominated by the central supermassive black hole, we have 

2(27r) 3 / 2 G 2 fi 2 n 



1v—\ r„ 



t max —4 
-Ul 



for U> < Qfcn 



(6) 



i.e, P x ,u is dominated by contributions from those objects near- 
est the central black hole and has a characteristic frequency 
dependence independent of the stellar distribution power law 
index. Thus, it is possible to make robust predictions for the 
scale of the jitter of Sgr A*, dependent primarily upon the max- 
imum density and typical masses of stellar-mass objects in Sgr 
A*'s immediate vicinity. 

Figure 1 shows the characteristic displacements, given by 
y/tjP XjU , for a handful of possible stellar/remnant core mod- 
els. In all but one the power spectrum follows it's asymptotic 
oc lu~ a form (corresponding to w~ 3//2 in the figure), despite the 
fact that n(r) is not a single power law and fi 2 (r) ranges from 
13.7M 2 at large radii to 1O 2 M 2 at small radii. What varies 
amongst the power spectra is primarily the normalization. 

The solid line shows our most optimistic case, with the 
stellar and black hole components given by Equations (3) 
and (4), respectively, truncated at the innermost stable circu- 
lar orbit (ISCO) of the central black hole, located at roughly 
1.2 x 10~ 6 pc. In this case the displacement spectrum is domi- 
nated by the black hole component, and produces fluctuations 
in excess of a /ias on 0.5 yr timescales. If, however, we ig- 
nore the remnant component altogether, the stellar component 
given by Equation (3) produces the dot-dashed line, implying 
that observable displacements will occur on 1 yr timescales. In 
both cases the typical displacements on year-long timescales 
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are comparable to those estimated in Reid & Bmnthaler (2004) 
(which find typical velocities of 0.3 km s" 1 and 0.07 kms" 1 , pro- 
ducing angular displacements of 8 /ias and 2 /j,as due to the rem- 
nant and stellar components, respectively 4 ). This presumes that 
we are justified in extrapolating the observed stellar density all 
the way down to the ISCO, over four orders of magnitude. If 
instead we force the stellar density to saturate at 0.01 pc, the 
smallest scale for which the stellar distribution has been mea- 
sured, we obtain the long-dashed line. Even in this case observ- 
able motions are present on sub-decade timescales. Finally, if 
we make the extreme pessimistic choice of truncating the stel- 
lar and black hole distributions at 0.01 pc we obtain the short- 
dashed spectrum, the curve of which is indicative of the lack 
of objects with sufficiently short periods to generate the short- 
timescale fluctuations. 

3.1.3. Searching for a Massive Binary Companion 

The hierarchical paradigm of galaxy formation also implies 
the presence of supermassive black hole binaries at the cen- 
ters of some galaxies. Initially the supermassive black holes 
sink toward the center of the combined galaxy as a result of 
dynamical friction due to stars in the reassembled bulge and 
drag imparted by gas funneled into the galactic center by the 
merger. However, this process becomes inefficient when the 
binary separation shrinks to a few pc, at which point the stellar 
encounters that carry away the orbital angular momentum be- 
come rare. This is before gravitation radiation is able to drive 
the system to merger within a Hubble time, resulting in the 
so-called "final parsec" problem (Begelman et al. 1980; Roos 
1981). 

A number of solutions to the final parsec problem have been 
suggested, including triaxial bulges (Merritt & Poon 2004), 
massive perturbers (Perets & Alexander 2008), and interac- 
tions with the merger enhanced gas density (Escala et al. 2004, 
2005). Nevertheless, in some fraction of systems long-lived su- 
permassive black hole binaries are expected (Begelman et al. 
1980; Yu & Tremaine 2003; Milosavljevic & Merritt 2003; 
Gualandris & Merritt 2009). For this reason searches for su- 
permassive binary companions have been performed in a va- 
riety of systems, primarily using spectroscopic methods (Bog- 
danovic et al. 2009; Boroson & Lauer 2009; Decarli et al. 2010; 
Liu et al. 2010). 

In the specific context of Sgr A*, massive companions and 
perturbers have been implicated for a very different reason : as a 
solution to the paradox of youth, i.e., the presence of the mas- 
sive, and therefore young, stars orbiting Sgr A* (Ghez et al. 
2003). Formation of the "S-stars" (typically B-type stars) in 
situ is thought to be precluded by the strong tidal forces from 
Sgr A*. This has lead to the development of models in which 
the stars are formed further out, where the restrictions are less 
severe, and then transported into Sgr A*'s vicinity (Gerhard 
2001). However, it is only possible to reach the small radii 
observed if the S-stars were members of a very tightly bound 
cluster, requiring a heretofore undetected massive component, 
the prime candidate therefore being an intermediate mass black 
hole (McMillan & Portegies Zwart 2003; Hansen & Milosavl- 
jevic 2003; Kim et al. 2004; Levin et al. 2005; Fujii et al. 2009; 
Merritt et al. 2009; Fujii et al. 2010). That is, the transport 
models imply the presence of a massive object located at a dis- 
tance comparable to that of the S-star orbits, roughly 10~ 2 pc, 

4 In Reid & Brunthaler (2004) the displacements due to remnants was twice 
as large as that reported here due to their assumption of twice as much mass in 
the black hole cusp. 
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Figure 2. Non-excluded region of the companion mass-orbital separation pa- 
rameter space that can be probed by astrometry with /tas precisions (cf. Hansen 
& Milosavljevic 2003). Thick solid lines show (from outside in) locations in 
the parameter space with 1 %, 5%, 25%, 67%, and 95% probabilities of not hav- 
ing been previously detected by previous astrometric studies at 7 mm with the 
VLBA (violated in the upper right) while exhibiting projected position shifts 
larger than a /ias at some point over the next decade (violated in the lower 
left). The dotted hatched area shows the approximate region described in the 
text that is ruled out by 7 mm VLBA studies. For reference, lines of constant 
maximum angular deviation are shown for circular orbits by the dashed lines, 
corresponding to 1 /ias (leftmost), 10/tas, and 100 /tas (rightmost). Gravita- 
tional merger timescales for circular orbits are shown by the thin solid lines, 
ranging from 10~ 3 yr (leftmost) to lGyr (rightmost). The period of the orbit 
(in the test particle limit) is shown along the top axis. Finally, the semi-major 
axes of the S-stars are indicated by the arrows. 

from Sgr A*. 

Thus, in addition to the stochastic buffeting from stellar-mass 
black holes and stars, Sgr A* may also move as a consequence 
of a massive binary companion. Whether or not this is observ- 
able astrometrically depends upon the orbital parameters, mass 
of the perturber and timescale of observations. Some set of 
these is already ruled out by the lack of an oscillatory signal 
in the 7mm proper motion studies (Reid et al. 1999b; Reid & 
Brunthaler 2004, 2005). If we assume circular orbits with pe- 
riods short in comparison to the observation time 5 , these rule 
out companion masses 

^ 8xl04 (l^)"' M - (?) 

(see also Hansen & Milosavljevic 2003). In the presence of 
eccentricity this expression is not a hard limit, since fortuitous 
orientations can result in small projected displacements. 

Adopting a simple two-body model (i.e., ignoring the many- 
body effects of the stellar/remnant cusp) we estimate the like- 
lihood of a given set of orbital parameters both satisfying the 
existing astrometric limits and executing deviations larger than 

5 For decade-long observations, this is true for a < 0.01 pc. Above 0.01 pc 
the orbital periods exceed 30 yr, and thus only a small fraction of the orbit 
would have been observed, decreasing the ability to probe these cases. Fur- 
thermore, at these distances the velocity is nearly constant, and thus degenerate 
with that due to the solar motion. 
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1 /ias over the next decade. This provides some idea of the rel- 
evant parameter space that /ias-astrometry can exclude in the 
near future. Figure 2 shows this likelihood, marginalized over 
orbital orientation and eccentricity (assuming a flat prior) as 
a function of orbital separation and companion mass. As ex- 
pected, large masses and orbital separations are already ruled 
out by Reid & Brunthaler (2005) (specifically, by requiring that 
the displacements be less than that observed perpendicular to 
the Galactic plane). While including eccentric orbits shift the 
location of this constraint, Equation (7) follows the contours 
of constant likelihood, falling nearly upon the 67% confidence 
level. At larger separations, the orbital period substantially ex- 
ceeds the timescale over which the 7 mm VLBA observations 
of Sgr A* have been performed, and their diagnostic ability 
decreases. At small masses and orbital separations the com- 
panion induces only sub-/ias shifts in the position of Sgr A*, 
and is thus not likely to be detectable by the EHT. Even when 
the companion mass is sufficiently large, highly eccentric or- 
bits pass within the ISCO, restricting the allowed parameter 
space further. Additionally, massive companions on small or- 
bits merge rapidly, though an intermediate mass black hole as- 
sociated with the S-stars would necessarily have entered the re- 
gion in the past Myr. At large separations the fluctuations due 
to stellar-mass black holes described in Section 3. 1 .2 provide a 
source of additional noise (not shown in Figure 2). However, 
this is confined to the extreme lower right of Figure 2 even in 
the most optimistic cusp models. Regardless, for reference, 
the dashed lines in Figure 2 show the lower limits achievable 
with 1 /tas, 10 //as, and 100 //as astrometric precision, assum- 
ing circular orbits. Thus, high-precision astrometry with the 
EHT should prove particularly effective for detecting massive 
companions to Sgr A*. 

3.1.4. Imaging Dynamical Accretion/Ejection Features 

The astrometric observations we have discussed thus far have 
focused upon detection of changes in the physical location of 
the supermassive black hole. However, astrometry is also sen- 
sitive to changes in the structure of the source. The observa- 
tion of radio, NIR and X-ray flares in Sgr A*, with charac- 
teristic timescales ranging from ~ 20 min to hours, is strong 
evidence for substantial variability near the black hole horizon 
itself (Genzel et al. 2003; Eckart et al. 2006; Meyer et al. 2006; 
Porquet et al. 2008; Marrone et al. 2008). Highly variable, tur- 
bulent, magnetized plasmas flows are a generic component of 
both accretion and jet models for Sgr A*'s emission, and thus it 
is not surprising that transient features arise near the black hole 
with dynamical times comparable to the orbital timescale at the 
ISCO. Sgr A*'s flares are therefore naturally explained within 
the context of transient structural changes to the emitting re- 
gion, e.g., orbiting hot spots generated by magnetic reconnec- 
tion events or shocks and confined within a disk, expanding 
blobs of energetic material buoyantly ejected from a disk, or 
hot spots launched outwards in a jet. If this is indeed the case, 
flares should produce observable temporary fluctuations in Sgr 
A*'s location. 

As new bright features emerge or move within the source the 
centroid of the image shifts, following the bright structures. 
The motion of the image centroid associated with orbiting hot 
spots, in the context of Sgr A*'s NIR flares, has been discussed 
at length (Broderick & Loeb 2005, 2006), and is the primary 
motivation for the GRAVITY instrument being developed for 
the Very Large Telescope Interferometer (see, e.g. Bartko et al. 
2009). At 7 mm the opacity of the accretion flow precludes the 
observation of spots orbiting near the ISCO. Nevertheless, it 



is possible to detect the motion associated with bright spots 
orbiting more than 3QGM/c 2 from the black hole using the 
VLBA (Reid et al. 2008). The EHT bridges the two regimes, 
and should access the innermost portions of the accretion flow 
with sufficient resolution to resolve any transient spots. 

Detecting structural changes does not necessarily require ab- 
solute astrometry. If the image is not dominated by the dynam- 
ical feature, approaches based upon closure quantities and/or 
the evolution of the source polarization are capable of detect- 
ing changes in the shape and orientation of the emission re- 
gion (Doeleman et al. 2009b; Fish et al. 2009). Alternatively, if 
the image uniformly brightens or expands, closure techniques 
can characterize the evolution of the emitting region. However, 
when a single transient feature dominates the emission, as has 
been observed to be the case for bright flares in the NIR, and 
this feature evolves structurally, absolute astrometry provides 
the only method to unambiguously detect the spatial variations. 

More importantly, since the dynamical timescale in Sgr A* is 
~ 10 min, the assumptions underlying Earth-rotation aperture 
synthesis are explicitly violated during flares. Specifically, we 
are no longer justified in assuming that the underlying source 
structure remains fixed throughout the night. For the EHT this 
is especially limiting; even when all existing and planned mm 
and sub-mm facilities are phased together, like the VLBA, the 
EHT will only be able to provide sparse instantaneous cover- 
age of the u-v plane. With a sufficiently well specified model 
for the variability in hand, it is possible to analyze the visibili- 
ties directly, eliminating the need to produce images altogether. 
However, in the absence of such a model phase referencing pro- 
vides the only means by which to relate visibilities at different 
times, producing approximate images of the time-averaged in- 
tensity. 

Furthermore, in the presence of a single dominant spot, 
phase referencing may be required to remove the fundamen- 
tal spot-disk degeneracy. That is, in the absence of an identify- 
ing feature that locates the black hole (e.g., the silhouette of the 
horizon, which may not be visible in the instantaneous images), 
fitting the self-calibrated visibilities alone will not be sufficient 
to distinguish a moving spot (either orbiting, infalling or out- 
flowing) about a fixed quiescent disk and an artificially fixed 
spot with an artificially moving quiescent disk. Of course, ab- 
solute astrometry naturally resolves this ambiguity. 

To address these issues, it is necessary to achieve sub-beam 
(< 30/ias) astrometric precision on the dynamical timescale 
of Sgr A* (~ 20 min). As long as a sufficiently strong phase 
reference is within ~ 1°, this should be possible with the EHT. 

3.2. Possible Phase References 

Phase reference candidates must satisfy a number of con- 
straints. Chief among these are brightness and compactness; 
the reference must be detectable in mm-VLBI observations. 
Most of what we describe in this section are efforts to char- 
acterize these properties. However, in addition the structure 
of the phase reference and Sgr A* must be stationary over the 
timescale of the experiment. In the case of the latter this is as- 
sisted by direct imaging, which can unambiguously identify the 
location of the black hole via it's silhouette against the nearby 
emitting material. For the phase references the case for stabil- 
ity is less clear, and ultimately, this will have to be empirically 
established. Here we will assume that this can be done. For 
quasars this is well justified, being located at extraordinary dis- 
tances. For masers this may not be, limiting their utility to 
short-term (days or weeks) uses only. Using multiple phase 
references provides a means to test this assumption and poten- 
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tially mitigate violations. 

3.2.1. Quasars 

The quasars J1745-282 (W56), J1748-291 (W109), and 
J1740-293 (GC 441) are located within 1° of Sgr A*, and have 
been used previously as phase references for proper motion 
studies of Sgr A* with the VLA (Backer & Sramek 1999). The 
first two have been successfully used in VLBI observations at 
wavelengths of 3.5cm and above (Reid et al. 1999b; Bower 
et al. 2001; Reid & Brunthaler 2004). We show the spectral en- 
ergy densities (SEDs) of these sources are in Figure 3. These 
are constructed from the literature and are therefore not con- 
temporaneous, varying over nearly 30 years (Isaacman et al. 
1980; Zoonematkermani et al. 1990; Bower et al. 1999, 2001, 
2002; Nord et al. 2004; Roy et al. 2005). As a consequence, 
the uncertainties reflect the larger of either that associated with 
the intrinsic measurement or that due to the variability implied 
by the multi-epoch observations. 

All of these are optically thin above 43 GHz. At frequencies 
above this transition we assume the SED is characterized by 
a power-law, which we extrapolate to obtain single-dish flux 
densities at 230 GHz and 345 GHz. Inferring the magnitudes 
of visibilities on long baselines then requires some assumption 
about the source structure; here we make the extreme simplify- 
ing assumption that the candidate references are point sources. 
Therefore, we compare the extrapolated flux densities directly 
to those presented in Table 1 . However, these may be consid- 
erable overestimates if substantial structure is not present on 
scales smaller than ~ 10 2 /ias in these objects. While the point 
source approximation may not be justified (e.g., the flat high- 
frequency spectrum of J1745-282 suggests the presence of a 
compact optically thick component, and thus extended struc- 
ture), the observed sizes of the candidate references at longer 
wavelengths are consistent with those of point sources, scatter 
broadened by the intervening interstellar electrons 6 . 

The steep spectrum of J1740-293, shown in the bottom panel 
of Figure 3, predicts fluxes well below the limits along any 
baseline at both 230 GHz and 345 GHz. This is unsurprising 
since J 1740-293 is already too dim at 43 GHz to be employed 
as a phase reference (Reid et al. 1999a). Thus, while we in- 
clude it here for completeness, it is almost certainly excluded 
as a potential phase references at shorter wavelengths as well. 
By comparison, both J1748-291 and J1745-282 are used as 
phase references for proper motion studies of Sgr A* (Reid 
et al. 1999a; Reid & Brunthaler 2004; Reid et al. 2008). 

J 1748 -291 transitions into the optically thin regime between 
1 .4 GHz and 8. 4 GHz, above which the SED appears well char- 
acterized by a power-law. As seen in the middle panel of Figure 
3, there is some uncertainty in the mm flux densities, deriving 
from the considerable uncertainty in the 22 GHz and 43 GHz 
fluxes. This arises in part due to the lower resolution of the 
VLA observations at 8 .4 GHz in comparison to the VLB A mea- 
surements at higher frequencies. While these agree within their 
respective uncertainties, including only VLBA data points pro- 
duces a flatter spectrum on average, and thus higher mm flux 
densities. Nevertheless, it is clear that J1748-291 is unlikely 
to be sufficiently bright to be used as a phase reference by the 
EHT-I. For the EHT-C, it is still too dim to be detected on the 

6 The interstellar electron scattering that broadens images of the Galactic 
Center is subdominant at mm wavelengths. Since the size of the scattering 
kernel, well described by an asymmetric gaussian, scales as A - , this is not 
true at longer wavelengths, where scattering necessarily induces correlated flux 
losses. 



Hawaii-CARMA baseline, though may achieve marginal de- 
tections on baselines including more than 10 ALMA antenna. 
For this reason we consider J1748-291 to be a marginal case. 

In contrast, J1745-282 is almost certainly bright enough to 
be used as a phase reference at 230 GHz for both the EHT-I 
and EHT-C, and 345 GHz in the case of the latter. The spec- 
tral index of J1745-282 is poorly constrained at present since 
the SED transitions to the optically thin regime near or above 
8.4GHz. More importantly, the 43 GHz flux estimates are lim- 
ited by atmospheric transmission effects, which bias the esti- 
mates towards lower values. As a result, the highest observed 
flux values may be better estimates of the high-frequency lumi- 
nosity of J1745-282 than the average value. Therefore, in the 
interest of completeness we show both, and their associated es- 
timates, in the left panel of Figure 3. In the former case the 
spectrum remains inverted at high frequencies, implying flux 
densities of ~ 0.1 Jy at mm wavelengths, well in excess of the 
various detection limits. In the latter case, the flux densities are 
more modest, ~ 30mJy, and comparable to the limits for the 
EHT-I. Nevertheless, both estimates imply that for the EHT-C 
J1745-282 is a strong mm- VLBI phase reference candidate. 

3.2.2. Masers 

Masers provide an attractive alternative due to their small in- 
trinsic size and large brightness temperature. Many molecular 
transitions exist that can produce masers at high frequencies, 
and those at mm and sub-mm wavelengths have only begun 
to be investigated (Humphreys 2007). Known sites of maser 
emission in the Galactic center include AGB stars (e.g., IRS 7, 
Menten et al. 1997; Reid et al. 2003, 2007) and nearby star for- 
mation regions (e.g., Sgr B2, Qin et al. 2008; Reid et al. 2009; 
Caswell et al. 2010). 

The SiO masers associated with AGB stars are believed to 
originate within 4-5 stellar radii (~ 8AU) of the star itself. 
However, typical AGB stellar radii are ~ 300/?©, implying an 
angular scale for the masing region of roughly 1 mas (Reid & 
Menten 2007). As a consequences, many of these sources will 
almost certainly be resolved out by mm- VLBI observations, 
eliminating them as candidate phase references. 

In contrast, H2O masers have been observed to have typi- 
cal extents comparable to an AU within star forming regions, 
corresponding to an angular scale ~ 10 2 /xas (Reid & Moran 
1981), and H2O maser spots an order of magnitude smaller 
have been observed surrounding late-type stars (Imai et al. 
1997). At 22GHz, H2O masers have already been used in 
VLBI astrometry observations within the Galactic center. Reid 
et al. (2009) report the detection of a number of strong maser 
sources in the Sgr B2 complex with typical sizes of ~ 0.3 mas. 
This is much smaller than the ~ 2 mas scatter-broadened size 
of a point source co-located with Sgr A* at 22 GHz, and thus 
it is unclear how much of the observed typical maser spot size 
arises due to interstellar electron scattering. As a consequence 
the Sgr B2 maser spot observations may be consistent with 
much smaller intrinsic sizes. High-frequency H2O maser tran- 
sitions at 325 GHz and 380GHz have been observed with fluxes 
similar to that of associated 22 GHz H2O masers surrounding 
late-type stars. 

As with quasars, the need to detect the masing region within 
the atmospheric coherence time produces a flux limit, below 
which the maser is not a viable phase reference candidate. Un- 
like continuum sources, however, the line-like nature of the 
maser prevents the EHT from leveraging large bandwidths to 
increase sensitivity. Rather, the signal-to-noise ratio is maxi- 
mized when the observing bandwidth matches the maser line 
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Figure 3. SEDs of radio-bright quasars within 1 ° of Sgr A* considered by Bower et al. (2001). The flux densities shown by the black circles were taken from the 
literature and are not contemporaneous (Isaacman et al. 1980; Zoonematkermani et al. 1990; Bower et al. 1999, 2001, 2002; Nord et al. 2004; Roy et al. 2005). Thus 
the error bars indicate the larger of the intrinsic measurement uncertainty or the variability amongst multiple data points at a given frequency. Green squares show the 
inferred flux densities at 230GHz and 345GHz. For J1745-282 the 43GHz flux measurement is highly variable, due presumably to observing conditions, and thus 
we also show the inferred high-frequency flux densities associated with the brightest 43 GHz measurement by the blue squares. For J1748-291 the blue squares show 
the mm flux densities inferred from VLBA observations alone. For reference the sensitivities of baselines consisting of stations with multiple telescopes at both ends 
are shown at 230 GHz and 345 GHz, by the open-orange and filled-red triangles (see the text for details). 



width, which for typical velocity widths of AV ~ 5kms _1 
are roughly 4MHz at 233 GHz. The corresponding detection 
threshold is then dependent upon the velocity width, with 

F " = F "°(5^)~" 2 < S > 

for which the coefficients F^ are given in Table 1 . These are 
typically on the order of 1 Jy for the EHT-I and 0.3 Jy for the 
EHT-C. By comparison, these are much smaller than the typ- 
ical flux of 40 Jy reported for the Sgr B2M spot observed at 
22 GHz by Reid et al. (2009). 

4. M87 

M87 is well known as a prototypical AGN, exhibiting a pow- 
erful relativistic jet. Recent analyses of its surrounding bulge 
have produced mass and distance estimates for the central su- 
permassive black hole of 6.6 ± 0.4 x 10 9 M Q and 17.9 Mpc, re- 
spectively (Gebhardt & Thomas 2009; Gebhardt et al. 2011). 
As a consequence, the angular size of M87's horizon is roughly 
70% that of Sgr A*'s, and thus M87 is an excellent second tar- 
get for the EHT 

4.1. Scientific Objectives 

The scientific objectives that can be addressed with short- 
wavelength astrometric observations in M87 are qualitatively 
different from those described for Sgr A* for two reasons. 
First, M87's much larger mass results in correspondingly sub- 
stantially larger dynamical scales. Unlike Sgr A*, M87 evolves 
slowly, with the period at the ISCO of > 5 days, and thus us- 
ing Earth-rotation aperture synthesis is well justified even for 
epochs of large variability. Similarly the time and mass scales 
for massive companions are increased by ~ 1.5 x 10 3 . Thus, 
decade long observations can only search for objects within 
~ 1.4 x l0 2 GM/c 2 , and more massive than ~ 10 7 M Q . How- 
ever, such a supermassive black hole binary will merge in less 
than 5Myr due to gravitational radiation alone, and is thus 
exceedingly unlikely. Finally, fluctuations in M87's position 
driven by stellar-mass black holes are essentially undetectable 



as a result of both, the increased timescale and the much larger 
mass of the central object. 

Second, in contrast to Sgr A*, M87 powers an ultrarelativis- 
tic jet. For this reason M87 has already been subject of a num- 
ber of high-resolution imaging efforts (Junor et al. 1999; Ly 
et al. 2004; Krichbaum et al. 2006; Kovalev et al. 2007; Ly 
et al. 2007; Walker et al. 2008). Nevertheless, the relation- 
ship between the observed jet emission and the black hole re- 
mains unclear (see, e.g. Marscher et al. 2008). By resolving 
the horizon, the EHT will unambiguously identify the central 
black hole. Registering the mm/sub-mm and radio VLBI im- 
ages will therefore allow the conclusive determination of the 
location of the black hole relative to the larger-scale radio fea- 
tures. Note that to achieve this, /ias precision is not required; 
the ability to position the radio VLBI and mm/sub-mm images 
is fundamentally limited by the astrometric precision at longer 
wavelengths, roughly 30 /ias. Positioning the black hole in this 
manner provides immediate insight into a variety of aspects of 
the relationship between the central supermassive black hole 
and the surrounding larger-scale emission, with associated im- 
plications for jet formation. In the remainder of this section we 
expand upon what can be learned by relating the mm and radio 
images. 

4. 1. 1. Locating the Counter Jet with Implications for Inclination and 

Acceleration 

Ly et al. (2007) and Kovalev et al. (2007) report tantaliz- 
ing evidence for the presence of a counter jet in M87. This 
arises from stacking multiple self-calibrated images, though it 
is probably not a calibration artifact. However, the interpre- 
tation of this feature as a counter jet hinges upon the identifi- 
cation of the core as the black hole location. In the absence 
of phase referenced VLBA observations of M87 showing sym- 
metric outflowing motions, this remains an assumption. 

The presence of a counter jet places severe constraints upon 
the inclination angle, O, and acceleration of M87's jet. Mak- 
ing the simplifying assumptions that the jet and counter jet are 
symmetric, with constant spectral indexes (Si, oc v~ a ), and only 
moving radially with Lorentz factor T, the observed brightness 
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ratio is approximately 

R _Ij„ ( 1+^cose \ a+1 „ ( 4 \ ~* 

Rb = 1^^{i-P C0S q) -{f=^J ' (9) 

where in the final expression the empirically well-justified ap- 
proximations r> 1 and 9< 1 were made. Note that this is 
a function of distance from the jet and in general will not be 
constant near the black hole. 

For large Lorentz factors, and thus at large distances, Rb, is 
a function of alone. The observed brightness ratio for the 
putative counter jet is ~ 15, and is roughly consistent with the 
inclinations of 30°^45° obtained by radio proper motion stud- 
ies (Reid et al. 1989; Ly et al. 2007). This is, however, much 
larger than the 20° upper limit implied by apparent motions 
of ~ 6c found by optical proper motion studies (Biretta et al. 
1999). There are many possible explanations for this discrep- 
ancy, including asymmetric jets, anisotropic emission, signifi- 
cant toroidal velocity components, significant jet structure and 
variable spectral indexes, as well as the misidentification of the 
counter jet. 

Near the black hole, within the jet acceleration region T can 
be considerably smaller than the asymptotic value at large dis- 
tances. This is consistent with radio observations of moving 
features, finding proper motions up to 2c within 2 mas from the 
core (Walker et al. 2008), in contrast to the fast moving optical 
features 500 times further away. When T < _1 , Rb becomes 
a rapidly varying function of distance from the core, providing 
an independent means to estimate the acceleration of the bulk 
flow within the jet as a function of radius. 

Utilizing the counter jet to study the jet acceleration and ori- 
entation is necessarily predicated upon an accurate determina- 
tion of the location of the central supermassive black hole for 
two reasons. First, robustly identifying a given image feature 
as being due to the counter jet is done most naturally by relating 
its position to that of the black hole and jet. Second, since the 
jet brightness evolves as it moves away from the launching re- 
gion, quantitatively comparing the jet-counter jet brightnesses 
even for the idealized, symmetric case, requires accurate dis- 
tances from the central black hole. Positioning the supermas- 
sive black hole within the core via EHT observations would 
remove any ambiguity regarding the black hole location, and 
would undergird the interpretation of the putative counter jet. 

4.1.2. Unambiguous Measurement of the Jet Opening Angle 

AGN jets are known to have remarkably small opening an- 
gles on large scales. This is set by both, the size of the jet 
launching region and the rapidity with which the jet is colli- 
mated. For M87, 7 mm images exhibit a radio core and a jet 
containing two limb-brightened features, with opening angle 
15°, commonly associated with the opening angle of the jet it- 
self (Junor et al. 1999; Ly et al. 2004; Krichbaum et al. 2006; 
Kovalev et al. 2007; Ly et al. 2007; Walker et al. 2008). How- 
ever, these limb-brightened features do not appear to converge 
upon the 7 mm core, intersecting roughly 2 mas east of it (Ly 
et al. 2007). One is then driven to one of three possible con- 
clusions, depending upon the location of the black hole and the 
breadth of the jet launching region. 

First, the 7 mm core does not indicate the position of the cen- 
tral black hole, and is instead a feature within the inner jet, e.g., 
as suggested by Marscher et al. (2008). In this case, it isn't 
clear why the region surrounding the black hole doesn't radi- 
ate. Nevertheless, the 7 mm jet morphology would be broadly 
consistent with jet formation theory. 



Second, the 7 mm core does indicate the position of the cen- 
tral black hole, and the jet rapidly collimates, begining from 
exceedingly large opening angles (> 60°) and reaching its 
asymptotic conical structure within ~ 0.3mas ~ lQ 2 GM/c 2 of 
the black hole (Junor et al. 1999). When the inclination exceeds 
the largest opening angle, projection of the jet on the plane of 
the sky significantly enlarges the apparent opening angle fur- 
ther, though this is insufficient to explain the rapid evolution in 
the opening angle by itself. In this case, M87 presents a seri- 
ous challenge to many jet formation theories, and jet formation 
simulations in particular. 

Third, the 7 mm core does indicate the position of the cen- 
tral black hole, and the 7 mm jet seen on mas scales is associ- 
ated with a different structure than the optical jet seen at large 
distances. This would naturally account for the systematically 
lower velocities observed in the radio images. In this case the 
transrelativistic unbound wind seen in many jet formation sim- 
ulations is a natural candidate for the structure responsible for 
the 7 mm jet feature. 

Locating the black hole with the EHT, and translating that 
location to the longer-wavelength radio images, conclusively 
addresses the distinction between the first two: either the 
black hole is located in the radio core or not. Furthermore, 
since the EHT promises to resolve the jet launching region, it 
will be capable of distinguishing between the compact ultra- 
relativistic jets and broad trans-relativistic winds (Broderick & 
Loeb 2009). By registering the mm-VLBI image to the longer- 
wavelength VLBA images, the relationship between the antic- 
ipated components of the outflow and the larger-scale features 
of the radio images may be determined. 

4.1.3. Identifying the Source of Jet Variability 

M87 is known to be variable on timescales of weeks, gen- 
erating bright knots which rapidly propagate away from the 
radio core. The apparent speeds of these knots vary consid- 
erably, ranging between 0.2c to 2.5c (Ly et al. 2007; Kovalev 
et al. 2007; Walker et al. 2008). In at least one case these have 
been accompanied by a TeV flare (Acciari et al. 2009; Wag- 
ner et al. 2009), implying that the source of the variability lies 
close to the black hole. However, where the variability origi- 
nates remains to be determined. Possible sources include varia- 
tions in the accretion rate, turbulence within the accretion flow, 
and instabilities within the jet itself (see, e.g., Giannios et al. 
2010). Monitoring M87 with coincident phase-referenced EHT 
and VLBA observations allows variable features to be followed 
from horizon to parsec scales. Those arising near the black hole 
will be clearly identifiable by their proximity to the silhouetted 
horizon. 

4.2. Possible Phase References 

Located well out of the Galactic plane, M87 does not suffer 
the interstellar electron scatter-broadening and extinction asso- 
ciated with studies of the Galactic center. This has facilitated 
the identification of a number of compact sources within it's 
vicinity. However, few of these are radio sources, and fewer 
still have been studied at VLBI resolutions. As a consequence, 
the primary uncertainty associated with the viability of poten- 
tial mm-VLBI phase references for M87 is their compactness. 
Nevertheless, we have identified a number of possible targets, 
which we list in order of decreasing viability. 

Located 1.52° away from M87, M84 (NGC 4374) is one of 
the few sources near M87 for which it is possible to find mul- 
tiple VLBI observations in the literature. As with many of the 
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Figure 4. SEDs of bright radio sources within 2° of M87, obtained from the literature. Left: SED of M84. The flux densities shown by the large filled circles 
correspond to observations with beam FWHMs smaller than 35". Those shown in blue are VLBI observations, with typical beam sizes of mas. The 23QGHz and 
345 GHz fluxes implied by all of these points and the VLBI points alone are shown by the green and blue squares (respectively). The small circles show all of the 
radio flux densities in the NED and VizieR databases associated with beams larger than 35", and are presumably indicative of the galactic non-thermal emission. 
Middle: SED of M89. Large light blue circles show VLBI observations, while large black circles show all other radio observations in the NED and VizieR databases. 
The green squares correspond to the 230 GHz and 345 GHz flux densities implied by all observations. For comparison, an upper limit at 347 GHz is shown by the 
arrow. Right: SED of M58. Large light blue circles show VLBI observations, while small open circles show all other radio observations in the NED and VizieR 
databases. The green squares are the 230GHz and 345 GHz flux densities implied by the VLBI measurements. In all panels, for reference the sensitivities of baselines 
consisting of stations with multiple telescopes at both ends are shown at 230 GHz and 345 GHz, by the open-orange and filled-red triangles (see the text for details). 



phase references we consider here, separating the contributions 
to the SED, shown in Figure 4, associated with the surrounding 
galaxy and the compact core is the primary difficulty. In the 
case of M84 we do this via the resolution dependence of the 
spectral index. The SED obtained from observations found via 
NED 7 and VizieR (Ochsenbein et al. 2000) with large beams 
(FWHMs > 35", small open circles, Heeschen & Wade 1964; 
Pilkington & Scott 1965; Pauliny-Toth et al. 1966; Gower et al. 
1967; Ekers 1969; Kellermann et al. 1969; Stull 1971; Witzel 
et al. 1971; Pauliny-Toth & Kellermann 1972; Slee & Higgins 
1973; Kellermann & Pauliny-Toth 1973; Condon & Jauncey 
1974; Heeschen & Conklin 1975; Sramek 1975; Wills 1975; 
Genzel et al. 1976; Veron 1977; Pauliny-Toth et al. 1978; Dres- 
sel & Condon 1978; Geldzahler & Witzel 1981; Large et al. 
1981; Kuehr et al. 1981; Bennett et al. 1986; Becker et al. 
1991; Gregory & Condon 1991; White & Becker 1992; Slee 
1995; Douglas et al. 1996; Condon et al. 2002; Vollmer et al. 
2004; Cohen et al. 2007) shows a characteristic steep power- 
law, almost certainly associated with the galactic synchrotron 
component. In contrast, the SED associated with observations 
with small beams (FWHMs < 35", large filled circles, Leeuw 
et al. 2000; Nagar et al. 2001, 2002; Leeuw et al. 2004; Ly et al. 
2004; Doi et al. 2005; Nagar et al. 2005), and VLBI observa- 
tions in particular (FWHMs ~ mas, large filled blue circles, Ly 
et al. 2004; Nagar et al. 2005), exhibit the flat spectrum indica- 
tive of AGN, extending from 5 GHz (the smallest frequency 
at which sufficiently small beam observations have been per- 
formed) to 670 GHz, above which dust emission begins to dom- 
inate (Leeuw et al. 2000, 2004; Doi et al. 2005). Below 43 GHz 
the small-beam and VLBI flux densities are similar, implying 
that in the radio, the small-beam observations are dominated by 
the emission on sub-mas scales. At far-infrared wavelengths, 
Doi et al. (2005) have excluded other potential sources of con- 
tamination (galactic synchrotron, free-free and dust emission), 
arguing that the dominant component is also due to the AGN it- 

7 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 



self. As a result, the existing 0.85 mm (350GHz) and 1.35 mm 
(221 GHz) flux measurements, obtained via the JCMT with 
sub-20" beams, provide a plausible estimate for the 230 GHz 
and 345 GHz fluxes that will be observed by the EHT; find- 
ing values in excess of 0.1 Jy (green-squares). Nevertheless, 
we also extrapolate flux estimates from the VLBI observations 
alone, finding estimates on the order of 0.07 Jy. In both cases 
these are considerably larger than the flux limits of the EHT-C 
at 230 GHz. For the EHT-I, M84 is still almost certainly suffi- 
ciently bright to be used as a phase reference at 230 GHz, and 
is marginal at 345 GHz. As a result, M84 is likely to be a strong 
mm- VLBI phase reference candidate. 

M89 is somewhat closer to M87, at an angular separation of 
roughly 1.19°. The radio SED of M89 is shown in Figure 4, 
in which all of the data points were taken from the literature 
(Heeschen & Conklin 1975; Sramek 1975; Dressel & Condon 
1978; Becker et al. 1991; Gregory & Condon 1991; Nagar et al. 
2001 ; Condon et al. 2002; Filho et al. 2004; Vollmer et al. 2004; 
Doi et al. 2005; Nagar et al. 2005; Healey et al. 2007; Braun 
et al. 2007) As with M84, there is a paucity of VLBI observa- 
tions, in this case only the one by Nagar et al. (2005), though 
the resulting core flux density is similar to VLA observations at 
the same frequency, differing by roughly 30% (cf. Becker et al. 
1991; Gregory & Condon 1991; Sramek 1975; Nagar et al. 
2001). This again suggests that the lower-resolution observa- 
tions are dominated by the compact core. Unlike M84, there is 
a substantial scatter in the high-frequency (~ 90 GHz) flux den- 
sities, possibly due to resolution, an interpretation supported 
by the fact that the relevant beam widths varied by an order of 
magnitude, with the smaller beam measuring the smaller flux 
density (e.g., Heeschen & Conklin 1975; Doi et al. 2005). This 
is unlikely to be responsible for the scatter at lower frequencies 
(near 10 GHz), where the higher-resolution observations pro- 
duce higher flux density estimates, though during a different 
epoch (e.g., Nagar et al. 2001; Vollmer et al. 2004). Therefore, 
to infer the 230 GHz and 345 GHz flux densities we fit a single 
power-law to all of the M89 radio data. This results in values 
that are consistent with the 347 GHz upper limit obtained by 
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Doi et al. (2005); sufficiently bright for the EHT-C, and only 
marginal for the EHT-I. Thus, we also consider this a strong 
candidate mm- VLB I phase reference candidate, though clearly 
less so than M84. 

The final candidate we show in Figure 4 is M58, again with 
the radio SED obtained from the literature (Heeschen & Wade 
1964; Sramek 1975; Dressel & Condon 1978; Gregory & Con- 
don 1991; Becker et al. 1991; White et al. 1997; Nagar et al. 
2001; Ulvestad & Ho 2001; Condon et al. 2002; Vollmer et al. 
2004; Anderson et al. 2004; Nagar et al. 2005; Dale et al. 2005; 
Doi et al. 2005; Gallimore et al. 2006; Krips et al. 2007). Sep- 
arated by 1.78° from M87, M58 is the further than M84 and 
M89, though still well within the isoplanatic angles at the rele- 
vant EHT stations on good nights. Unlike M84 and M89, suf- 
ficient VLB1 observations exist for M58 to completely charac- 
terize its radio spectrum, from 1.7 GHz through 43 GHz (An- 
derson et al. 2004; Nagar et al. 2005; Krips et al. 2007). These 
allow a reasonably precise estimate of the mm-VLBI flux den- 
sities, giving values comparable to the sensitivity limits of the 
EHT-C at 230 GHz, and well below them along the Hawaii- 
CARMA baseline at 345 GHz. As a consequence, this is a 
marginal phase reference candidate. 

There are a number of additional bright sources near M87, 
though none are well studied at sufficiently high resolutions 
to make definitive statements regarding their millimeter fluxes 
on sub-mas scales. Within 1°, the most promising candidates 
are J1229+114 (PKS 1227+119, VPCX 27) and J1230+114 
(VPCX 28). Extrapolating from the existing radio data (Shim- 
mins et al. 1975; Large et al. 1981; Wright & Otrupcek 
1990; Becker et al. 1991; White & Becker 1992; Ledlow & 
Owen 1995; Gregory et al. 1996; White et al. 1997; Laurent- 
Muehleisen et al. 1997; Condon et al. 1998; Reich et al. 2000; 
Cohen et al. 2007) gives estimates for the millimeter fluxes 
comparable to, or in excess of, the detection limits of the EHT- 
C, though in most cases comparable to the limits for the EHT-I. 
However, these conclusions are contingent upon a substantial 
fraction of the radio flux arising from sub-mas scales. There 
is some evidence for the radio emission in J1229+114 being 
primarily associated with the AGN, as opposed to galactic syn- 
chrotron; J1229+1 14 has been classified as an FRI based upon 
the 1.4GHz morphology as measured at 5" resolutions by the 
FIRST survey, performed with the Very Large Array (White 
et al. 1997; Gendre & Wall 2008). However, at 5" it is already 
apparent that the source is extended, with the unresolved core 
constituting approximately 10% of the overall flux (White et al. 
1997). Similarly, Flesch (2010) assigns a 41% probability that 
J1230+1 14 is associated with a quasar based upon radio-X-ray 
associations. Furthermore, it is the only radio source within 
5' of MC2 1227+120, classified as a quasar in Veron-Cetty & 
Veron (1983), despite the lack of this source in the Molongo 
catalogs (cf. Large et al. 1981). Nevertheless, since neither ob- 
ject has been observed with resolutions exceeding 5", whether 
or not J1229+1 14 and J1230+1 14 are suitable phase references 
remains highly uncertain. 

Finally, many of the science objectives described in Section 
4.1 require only registering the mm-VLBI and cm-VLBI refer- 
ence frames. As a consequence, many of these can be achieved 
with astrometric accuracies of ~ 0.1 mas. Furthermore, due to 
the long dynamical timescale of M87, it is possible to com- 
bine observations over many nights to reduce the astromet- 
ric error. However, as shown in Appendix A, the constraint 
upon the proximity of the phase reference is dependent upon 
the astrometric precision desired and the length of the obser- 
vation. Specifically, Equation A12 implies that 0.1 mas pre- 



cisions may be achieved at all relevant EHT stations over the 
course of a single night, corresponding to hundreds of indepen- 
dent 10 s measurements, for S < 20°. While this is sufficiently 
large that it violates some of the assumptions made in deriv- 
ing Equation A 12, it is nevertheless clear that angular separa- 
tions much larger than a few degrees are permissible for this 
purpose. Thus, we include a fifth possible reference, the well- 
studied bright quasar 3C 273, located 10.3° away from M87. A 
sufficient number of VLBI flux measurements exist for 3C 273 
to produce an SED from these alone (Krichbaum et al. 1997; 
Lonsdale et al. 1998; Kellermann et al. 2004; Zavala & Tay- 
lor 2004; Horiuchi et al. 2004; Lister & Homan 2005; Kovalev 
et al. 2005; Cara & Lister 2008; Lee et al. 2008). They include 
the observation by Krichbaum et al. (1997) at 215 GHz which 
found a flux density of 9.2±0.6Jy. As a result, there is no 
question that 3C 273 is bright enough at millimeter and sub- 
millimeter wavelengths to be used as an astrometric reference. 

5. CONCLUSIONS 

The EHT promises to usher in an era of /ias black hole as- 
tronomy. This takes the form first of ~ 20 /ias-resolution im- 
ages, allowing careful study of the morphology and dynamics 
of accretion flows and jet formation on sub-horizon scales for 
a handful of supermassive black holes, including Sgr A* and 
M87. Because these images resolve the horizon, they unam- 
biguously locate the black hole via its silhouette against nearby 
emission. However, it also allows high-precision astromet- 
ric observations, relating the positions of Sgr A* and M87 to 
nearby, bright, compact sources. This is only possible because 
many of the EHT VLBI stations are arrays, allowing coincident 
observations of the source and a reference without having to re- 
point the antennas on the atmospheric coherence time, roughly 
10 s. On a typical night, atmospheric turbulence limits the pre- 
cision to 4 /ias over 10 min and 1 /ias over many hours. 

Necessary for /ias-precision astrometry with the EHT is 
nearby, bright point sources. This is complicated by both the 
short wavelength (few objects are bright in the mm) and the 
very high angular resolution of the EHT (many objects are re- 
solved out). Nevertheless, with the anticipated 4GHz band- 
widths a variety of sufficiently bright possible phase references 
exist near Sgr A* and M87. While in most cases we cannot 
determine the mm-size of these objects, in all cases they have 
not been resolved out at 7 mm. Note that we have concerned 
ourselves solely with objects that can be detected individually 
on timescales shorter than the atmospheric coherence time, and 
thus the sensitivity /brightness limits can be relaxed if these ob- 
jects are phase referenced by the primary source (e.g., Sgr A* 
or M87). 

The kinds of science that /ias-precision astrometry enables 
depends critically upon the object under consideration. Many 
of the processes that potentially affect the position of the cen- 
tral, supermassive black hole depend upon its mass, while the 
questions of interest depend upon its environment. In all cases, 
however, the unique ability of mm-VLBI to conclusively locate 
the position of the supermassive black hole is used. 

Sgr A* is constantly buffeted by objects within the evolved 
stellar cusp. Even with extremely pessimistic assumptions 
regarding the content and structure of the cusp, over year 
timescales, this produces /ias displacements in Sgr A*'s po- 
sition. While the presence of these displacements is robust, the 
magnitude is sensitive to the properties and density of the ob- 
jects within the cusp, and thus provides a means to probe these. 
Detection of the expected stellar-mass black hole cusp would 
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have far reaching implications for the event rates of existing 
and future gravitational wave searches. 

Sgr A* is also potentially perturbed by massive binary com- 
panions. Frequently invoked to explain the presence of young, 
massive stars within 10 2 AU of Sgr A*, an intermediate black 
hole orbiting Sgr A* has an infall timescale (dominated by 
gravitational radiation) that is long in comparison to the life- 
times of the young stars. If such a binary exists it will induce 
periodic displacements in Sgr A*'s position, the size of which 
depend upon the object's mass and orbital parameters. While 
some portion of the available parameter space has previously 
been excluded by 7 mm astrometric studies, the /ias-precision 
astrometry made possible with the EHT is capable of elimi- 
nating the vast majority of the remaining parameters with sub- 
decade observations. 

Efforts to image Sgr A* must necessarily contend with its 
sub-hour, large-amplitude variability. This formally precludes 
Earth aperture synthesis, which assumes that the source is fixed 
throughout the night and is not amenable to self-calibration 
techniques, which smear dynamic images in an ill-defined fash- 
ion. Phase reference observations provide a way in which to 
faithfully relate visibilities during variable episodes, providing 
a time-averaged image of the Sgr A*. In the absence of an a 
priori accurate model for Sgr A*'s flares, this makes it possible 
to identify the qualitative features of the relevant processes. 

Due to its much larger mass, what can be probed via astrome- 
try forM87 is very different than what is accessible for Sgr A*. 
In M87's case, perturbations induced by the stellar/remnant 
cusp are negligible, as are displacements due to potential binary 
companions. Instead, the utility of astrometry with the EHT is 
the ability to relate the mm-VLBI image, in which the black 
hole is clearly identifiable due to its silhouette against nearby 
emission, and the large-scale structures in longer-wavelength 
images. This provides a way in which to conclusively locate 
the black hole, eliminating ambiguity in the identification of a 
putative counter jet, inclination, opening angle, nature of the 
radio core, and relationship between the black hole and the ra- 
dio jet. Moreover, during dynamical periods, with the EHT 
it is possible to follow the creation and propagation of radio 
knots from sub-horizon to parsec scales, shedding light upon 
the source of the observed variability in M87 specifically, and 
AGN more generally. 
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APPENDIX 
A. MILLIMETER ISOPLANATIC ANGLE 

Understanding atmospheric propagation effects is critical to 
high-fidelity image reconstruction at high frequencies. As 
a result, considerable effort has been expended theoretically 
and empirically characterizing the phase fluctuations associ- 
ated with spatially varying indices of refraction in the turbulent 
atmosphere. A general treatment can be found in Thompson 
et al. (2001); here we collect some results that are of interest to 
mm-VLBI astrometry. 



We wish to estimate the limits imposed by atmospheric 
propagation effects upon the accuracy with which the rela- 
tive locations of nearby sources may be determined at 230 GHz 
(1.3mm) and 345GHz (0.87mm). To do this we will make a 
number of simplifying approximations and convenient assump- 
tions regarding the structure of the intervening atmosphere and 
the sources of interest. We begin by defining a model for the 
atmospheric fluctuations, and the associated phase errors along 
nearby lines of sight. After doing so, we will specify an ideal- 
ized astrometric experiment, and compute the associated errors 
induced by the forgoing atmospheric model. 

The atmospheric model is defined by fluctuations on top of 
some mean structure. We make the common assumption of a 
plane parallel atmosphere, imparting an average phase delay. 
Superimposed upon this are gaussian phase fluctuations, as- 
sumed to be due to a thin layer within the troposphere (fixed at 
an altitude of roughly L ~ 1 km). The resulting phase pertur- 
bations are then characterized by the mean value and the spa- 
tial correlation function, commonly referred to as the structure 
function. These are defined by their zenith values: 



$ z = (<fc(x,f)>, 
P (|d|) = ([<Mx+d,f)-<fc(x,f)] 2 
6l(x,f))-$I=^(oc) 



(Al) 



^2 



where we have given the auto-correlation function a unique 
symbol despite being degenerate with the V^d), and d will 
correspond to interferometer baselines. Each of these quanti- 
ties is a property of the particular site and may be empirically 
determined. For example, $z may be obtained from various 
atmospheric models (e.g., Niell 1996; Boehm et al. 2006a,b), 
with typical values of (600cm/A)rad (though this varies by 
roughly a factor of two between sea level and 5000 m, the al- 
titude of the ALMA and Hawaii sites). At small spatial sepa- 
rations the structure function has been found to be well fit by 

a power law, V^id) ~ (2ir<7do /A) (d/ 1 km) 2 ' 3 rad with values 
for ovo and f3 given in Table 2 for the relevant locations. At 
large d, V^id) asymptotes to a fixed value, 2E 2 4 . Typical val- 
ues for are 0.01<I>z — (6cm/A)rad, is dominated by large- 
scale variations in the wet component and varies substantially 
on hourly timescales. 

The phase perturbation along a general oblique line of sight, 
defined by a position on the sky a, associated with zenith angle 
zx.a at location x, in a plane-parallel geometry is given by 

(j> a (x, a, t ) = <fc(x+Ltanz x ,a , sec z x ,a • (A2) 
From this we can construct an angular structure function: 

X> (a,£)= ([(j) a (x,a! + 5,t)-<f) a (x,a,t)] 2 ) 
- <I?| sec 2 z tan 2 z 5z 2 



X> (c/)sec 2 z 



1 + tanz<5z + 



1 + sin z 2 
2cosz 



(A3) 



where we have suppressed the subscripts on z and Sz, which 
should be understood to be uniquely defined in terms of x, a, 
and 6, 



+ E 2 sec 2 ztan 2 z<5z 2 + 0(fe 



d = L\l [tan(z+cfc)-tanz] 2 + \S\ 2 -Sz 2 



(A4) 



LA/(sec 4 z-l)& / +|<5| 
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Table 2 

Phase Error Structure Function Parameters and Isoplanatic Angles at Potential EHT Stations 



Site 


Telescopes 


Alt. 




6f° 


GHz 


5 345GHz 


Reference 






(m) 


(mm) 






Mauna Kea c 


JCMT, SMA, CSO 


4070 


0.4-2.7 0.75 


0.6°- 


-7.7° 


0.4°^1.5° 


Masson (1994) 


Chile c 


ASTE, APEX, ALMA 


5000 


0.3-1.5 


1.6° 


-11° 


1.0° -6.6° 


NRAO (1998) 


Plateau de Bure c 


IRAM Interferometer 


2552 


0.3-0.7 1.1-1.9 


5.3°- 


■9.6° 


4.3°-6.7° 


Olmi & Downes (1992) 


Cedar Flat d 


CARMA 


2196 


0.8-3.8 


0.5°- 


-3.3° 


0.3°-2.0° 


Woody et al. (2004) 



Note. — Adapted from Thompson et al. (2001), Table 13.3. 
a Defined at a baseline of 1 km. 

b Where no value was given, /3 = 5/6, consistent with Kolmogorov turbulence, was assumed. 
c Ranges indicate good/bad observing conditions (e.g., nighttime/daytime). 
d Ranges indicate before and after water vapor radiometry has been employed. 



and we have made use of the Taylor expansion of sec(z+ 5z) to 
2nd order: 



sec(z+&) = secz 



1 +tanzfo+ 



1 +sin 

2cosz 



'■■Sz 2 + o(Sz 3 



(A5) 

Unlike T>,p{d), which is completely defined by the small-scale 
structure of the troposphere, V^ia.^) has a contribution due 
to large-scale atmospheric variations as a result of the differ- 
ent path lengths along the two lines of sight, resulting from 
the slightly different zenith angles. Which term dominates in 
Equation (A3) depends upon the zenith angle, how strongly 
nearby atmospheric patches are correlated and the magnitude 
of the large-scale atmospheric fluctuations about the average 
phase delay. 

To estimate the limits upon the accuracy this implies we 
consider a simplified astrometric experiment. Specifically, we 
will assume that the two sources of interest are monochromatic 
point sources separated by more than a beam width, though 
extension to more complicated source structures is straightfor- 
ward. Furthermore, we focus only upon atmospheric phase 
errors, and explicitly neglect all other contributions. Finally, 
we will assume the experiment consists of two sets of antenna 
pairs, one pair located at x = and the other pair at some po- 
sition x = u, sufficiently far away such that the atmospheric er- 
rors are completely uncorrelated. With these assumptions, the 
voltages induced by the radio waves emitted by the two point 
sources, located at positions a and a + S, at a pair of neigh- 
boring antennas at some position x, each viewing one of the 
sources, are 

ojt-\ — — • a. + <p a (x,a,t) 
A 

u)t + —j— ■ (a + S) + <p a (x,a+ S,t) 

(A6) 

where Eiji are two complex constants, <f> a (x,ot,t) are the ran- 
dom phase perturbations imparted by propagation through the 
atmosphere defined above, and the remaining term is the geo- 
metric phase delay that we seek to measure. Cross-correlating 
the antenna voltages along the long baselines produces corre- 
lation peaks for each source at time delays of 



V , /4 (x,0 = £'/exp|-/ 
V B (x,t) = E n exp\-i 



2 ™ J. , s 



T B = 



A 

27TU 



(A7) 



(a+S) + tf> a n(a+S,t), 



sumed that the atmospheric phase delay varies slowly in com- 
parison to the timescales over which the visibilities are mea- 
sured. From these the relative positions of the point sources 
may be extracted in the normal way, obtaining an estimate 
of the projection of the angular distance between two sources 
along u of 



u-<5 



A 

2-7TM 
:fi'5 + 



(t b -t a )-A 

A 
2™ 
A 



[(/tenia + S,t)~ 4>a,Vi{OL, t )] - A 



(A8) 



= u-S+ - — {[0 a (u, a+ S,t) - (j) a (n, a, t )] 
liru 

-[<t> a (0,a + d,t)-(f> a (0,a,t)]}-A, 

where A = <3> z secz tan z<5z| u - <E>z secz tan z<5z|q is the angular 
shift due to the differential phase delays at the two stations 
associated with the mean atmospheric model. The u ■ d b s is 
a gaussian random variable centered upon the true value and 
with a variance due to the atmospheric phase fluctuations of 



4 



(A9) 



This has contributions from both, the small-scale, rapidly vary- 
ing structures that limit the atmospheric coherence time (the 
V^(d) terms in Equation (A3)), and large-scale, slowly evolv- 
ing fluctuations (responsible for the terms in Equation 
(A3)). Because f3 is typically close to unity, which dominates 
the astrometric uncertainty is only weakly dependent upon the 
angular separation of the two sources. In the remainder of this 
section we discuss the two kinds of terms separately. 

The error term associated with the rapidly varying, small- 
scale structures due to a single station is given by 



^ A a d Q 
^ secz 

u A V 1km 



sec z 5z 



(A10) 



where <fr a . u (a,t) = <fi a {u, a,t)-<j> a {0,a.,t), in which we have as- 



where we have chosen 5z = \S\, resulting in the maximum a s for 
a given angular separation. For a values of ado, and z typical 
of the stations likely to be employed by the EHT, at 1.3 mm, 
a s < 0.5(A/t/)(L/lkm)' 3 ( ( 5z/l )' 3 . 

Since a s is dominated by contributions from the most rapidly 
varying structures in the atmosphere, observations separated in 
time by more than the atmospheric coherence time probe inde- 
pendent realizations of the atmospheric turbulence. As a con- 
sequence, during a few hour observing period it is possible to 
reduce a s by an order of magnitude. Thus, the limit upon 5 im- 
posed by a given astrometric tolerance, a a, and independent 
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observations is 



<5<£ = 57°cos 2 z 



lkm 



Ncosz 



UOA_ 



(All) 



Limits for N = \ and a& = X/u are listed in Table 2 for 1 .3 mm 
and 0.87 mm for z = 70°, typical of the zenith angle for Sgr 
A* as seen from North American sites (the most pessimistic 
scenario we consider). In terms of these, we have 



2.9cosz 



, 0.12sec-z — 
N \ OqJ u 



A 



(A12) 



The error associated with the large-scale, slowly evolving 
component at a single station is given by 



cri 



X 

U 27T 



sec z tan z Sz. 



(A13) 



For Sgr A* viewed from North America, typically 07 < 
(X/u)(5z/l°), can be substantially smaller at high-altitude 
sites, and is a strong function of zenith angle. 

Unlike the small-scale component, the long timescale over 
which the large-scale structures vary means that only a hand- 
ful of independent realizations will be encountered on a given 
night. Thus, substantially reducing 07 by integrating over long 
times requires observations spanning many days at least. How- 
ever, as we discuss in the main text, the large spatial and tem- 
poral scales of these fluctuations make possible a number of 
schemes to measure the additional path delay directly and sig- 
nificantly reduce the unmodeled S^,. For this reason the un- 
certainty in the angular separation due to unmodeled slowly 
varying and rapidly varying components, 



2 21 2 

a 8 = a s lo + °l 



+ of .,+ 07 



(A14) 



are likely to be similar. 



B. MULTI-REFERENCE ASTROMETRIC CALIBRATION 

Here we describe explicitly a scheme to measure and re- 
move the contribution to the astrometric uncertainty arising 
from slowly varying, large-scale atmospheric fluctuations (i.e., 
those responsible for the 07 terms in the preceding appendix). 
The method utilizes multiple nearby reference sources to mea- 
sure the slowly-varying phase fluctuations, which may then be 
explicitly accounted for in the astrometric measurement of in- 
terest. The viability of the scheme is predicated upon the fol- 
lowing assumptions: 

• The existence of a set of at least three reference sources 
near the target source with different zenith angles at each 
VLBI station. Note that the limitations upon the separa- 
tion derived in Section A may be substantially relaxed 
since both 0/ and a s are both roughly proportional to \8\, 
and thus the limit upon the fractional accuracy imposed 
by rapidly-varying, short timescale atmospheric fluctua- 
tions is only weakly dependent upon separation. In prac- 
tice |<$| is limited by the size of the region over which 
the large-scale phase perturbations are uniform, which is 
typically many degrees. 

• Over the portion of the sky of interest, the large-scale 
atmospheric phase delays at each site is well modeled by 
a single-parameter model (e.g., the secz model we will 
employ here). 



• The reference sources do not evolve intrinsically or rela- 
tive to each other. More specifically, the angular separa- 
tion between each of the reference sources is fixed. 

• The "true" angular separations of the reference sources 
can be measured with the desired fractional accuracy. 
This may be accomplished, e.g., by careful atmospheric 
modeling and long-time integrations. 

To illustrate the procedure, and estimate the resulting calibra- 
tion uncertainties we consider observations of three reference 
sources on a single, static baseline. We use the letters A, B, and 
C to indicate a particular source, and numbers 1 and 2 to indi- 
cate the locations of individual mm- VLBI stations. For short- 
hand we define A secz| B1 = secz^i -seczsi to be the difference 
between the secant of the zenith angle of reference sources A 
and B as seen from station 1, with similar expressions for all 
sources and both stations. In terms of these, following the pre- 
vious appendix, from a single observation we have an estimate 
of the separation of reference sources A and B of: 



U12 • OAB.obs = U12 • $ab + AS s j.B 
+ — ^A$, Asecz|£- — 



■A$2 Aseczl^ , (Bl) 



where we have separated the contributions to the phase fluc- 
tuation into modeled large-scale and the remaining unmodeled 
small-scale fluctuations (A8 S ). The A$i.2 are the unknown 
zenith phase delays implied at stations 1 and 2, respectively, as- 
sociated with the large-scale atmospheric perturbations. Since 
6ab is known, or may be measured by comparing data sets 
taken over many days, we may invert this to obtain an expres- 
sion for the unknown A<I>i .2: 



A$i Asecz|o!-A$ 2 A seed™ 



= ^r- ki2 • (8. 



X 



MB 



-AS 



sAB\ 



A$! Asecz|*i-A$ 2 Asecz^ 

2n r , , 1 (B2) 

= — |U12 • (OflC.obs - 8 B c) ~ AS Si BC\ 



A$j Asecz|^|-A$ 2 Asecz|" 

= — [Ul2 ' (8cA.obs- 8 C a) ~ AS SiC a] , 

where we have added analogous terms for the other source 
pairs. In principle, these are degenerate, providing only two 
linearly independent equations for A$i and A<£>2- In practice, 
since the AS S terms are dominated by the small-scale atmo- 
spheric structure, this degeneracy is weakly broken, providing 
an external estimate of the magnitude of the A5 S terms. Con- 
sidering only the first two, we obtain an estimate for A$i of 

A$i = 

27TUi2 Aseczl^ (<W>bs-<*Afl) - Aseczl^ (foc.obs - osc) 



Asecz|^2 Asecz|*J- Asecz|^[ Asecz|*2 



(B3) 



with a similar expression for A$2 obtained by transposing 1 
and 2 everywhere. Comparing this to Equation (B2), assuming 
that the errors in the 8ab.bc ma y be neglected in comparison to 
AS S (the errors in the zenith angle estimates result in higher- 
order corrections, and thus always may be ignored), we find a 
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variance in A$i induced by the small-scale, rapid variability 
of 



2 



(Asecz|^) o 1 sAB +( K 



w +\ Aseczl^ 



's.BC 



Asecz| B2 Asecz| cl ■ 



Asecz|^[ Asecz|*2 



2 ' 



(B4) 

with a similar expression for CTa*, ■ Note that if all of the 
reference sources lie at similar zenith angles, and have sim- 
ilar separations in zenith angle from each other, this reduces 
to <ta$ 12 ~ &s> and thus the fractional precision of the large- 
scale phase corrections is comparable to that associated with 
the small-scale, rapidly varying fluctuations over hour-long in- 
tegrations, roughly 1 /ias. The accuracy will be limited by how 
well the true separations of the reference sources can be mea- 
sured, presumably by averaging over long times, and the valid- 
ity of the assumptions underlying the scheme. 

C. POSITION JITTER POWER SPECTRUM 

Supermassive black holes settle via dynamical friction into 
the center of their galactic potential well on timescales short in 
comparison to the Hubble time (Chatterjee et al. 2002). Never- 
theless, fluctuations in the surrounding stellar/remnant 8 distri- 
bution induce random displacements of the central black hole 
similar to Brownian motion. Here we derive an estimate of 
the displacement power spectrum. For simplicity we assume 
that the surrounding distribution of gravitating bodies is spher- 
ically symmetric, with number density n{r) and mass function 
(f)(m,r), all orbits are circular, and that dynamical collective 
N-body effects are unimportant, i.e., the fluctuations about the 
mean density may be approximated by Poisson statistics. This 
last assumption may not be strictly justified in the spherical 
case, though is in the more generic triaxial situation. 

In the vicinity of the supermassive black hole, we may ap- 
proximate the black hole's motion by 



x + 7X+WQX = a(f), 



(CI) 



where x is the location of the black hole, a is the fluctuating ac- 
celeration induced by perturbations in the surrounding stellar 
distribution, uiq = d 2 $/dr 2 is the natural frequency of oscilla- 
tion in the bottom of the galactic potential well, and 7 is the 
damping associated with dynamical friction. We may estimate 
wo in terms of the stellar density at the galactic center: 



4tt 



(C2) 



where po is the average stellar density within the region of in- 
fluence. Thus, cjo is a proxy for the stellar mass density within 
the core. The dynamical damping rate is approximately given 
by 

7 ~ 4v / 2^1n AG 2 M^ (C3) 
CT o 



where In A is the Coulomb integral, M is the mass of the central 
supermassive black hole, and <tq is the stellar velocity disper- 
sion in the stellar core (Chatterjee et al. 2002). Typically, we 
will be interested in the case when to ^ ljq, 7, though we have 
kept all terms in the interest of generality. 

8 Henceforth, for simplicity we will refer to the "stellar" distribution, "stel- 
lar" mass, and "stars" though this should be understood to include evolved 
objects like stellar mass black holes as well. 



Since the system is spherically symmetric, we will focus 
only upon the displacement along the x direction. In this case, 
the power spectrum of these displacements is then given by 



Px,uj — (l-^u 



(ui 2 - Wq) 2 + 7 2 w 2 



(C4) 



where P fl w is the power spectrum of the fluctuating accelera- 
tions. Thus, computing the spectrum of displacements is re- 
duced to computing the spectrum of the fluctuating forces. 

In Chatterjee et al. (2002) it was assumed that the fluctua- 
tions were completely uncorrected in time, i.e., that P a ^ was 
constant. For virialized systems this may be derived directly 
via the fluctuation dissipation theorem. However, this is likely 
to be a poor approximation in the case of interest, where we 
are interested exclusively in short-timescale variations in the 
position of the supermassive black hole. Thus we approach 
the problem phenomenologically: given a density profile we 
estimate the P a M associated with Poisson fluctuations in the 
number density of stars. 

We begin by considering a thin shell of width Ar centered 
upon a radius r. Within this shell there aie N = 47rr 2 «(r)Ar 
objects, with mass function <fi(m, r), randomly distributed. The 
induced acceleration along the x-axis is then 



EGnijXi 
r 3 



(C5) 



from which it is obvious that (a) = 0. The variance in the ac- 
celeration is 

N n 2 n 2 N 

1 2\ u niimjXiXj u \ ■> . . . 



i,j=i ij=\ 

j=i 



(C6) 



where (xjXj) = r 2 Sij/3 was used and /i 2 (r) = dmm 2 (j)(m,r). 
Over densities within the shell will persist until the orbits of 
stars within the shell can re-randomize the stars. Therefore, 
unlike Chatterjee et al. (2002), we expect these perturbations to 
be correlated for a time comparable to 1 /f2t(r), where fifc(r) is 
the Keplerian angular velocity at r. Thus, the auto-correlation 
function of the acceleration fluctuations is approximately 



Rf T T = (a(t)a(t+T)) ~ (a 2 )e-^ T2 ' 2 = 9h^ Ne ~^T 2 l2 (C7) 



This implies a power spectrum associated with the shell of stars 
of 



lArr 2(2^GV » :V/m ; Ar 



r 2 Q k 



(C8) 



Since we have already assumed that the positions of individ- 
ual stars are uncorrected (when we assumed Poisson statistics 
within a shell), the power spectrum of the entire distribution of 
stars may then be constructed by summing multiple shells, i.e., 



P a .u= / dr 



2(2tt) 3 / 2 GV n 



/ml 



(C9) 



which may be inserted into Equation (C4) to obtain the desired 
power spectrum. 
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To make further progress we must choose <j)(m, r) and n(r). 
From these we may then construct the Sl^: 

- GM G f°° f°° 
\l k (r) = — — + — / dr I amm(p(m,r)n(r) , (CIO) 
r r Jo Jo 

and thus P x u . 
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